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We investigate non-minimal F'^-type couplings of electromagnetic fields to gravity. We derive 
qH' the field equations by a first order variational principle using the method of Lagrange multipliers, 

^jj^i Then we present various static, spherically symmetric solutions describing the exterior fields in the 

vicinity of electrically charged massive bodies. 
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I. INTRODUCTION 



Non-minimal couplings between electromagnetic fields and gravity can be induced in regions 
of space-time such as the neighborhood of black holes or neutron stars where high intensities 
of electromagnetic and/or gravitational fields are present. In such cases, gravitational fields 
behave as a non-linear medium in which the electromagnetic fields propagate. On the other 
hand, very intense electromagnetic fields may cause new gravitational effects. 

Non-minimal couplings of the form RF^ were investigated to obtain more information on 
the relationship between space-time curvature and electric charge conservation 

• They 

were also derived by Kaluza-Klein reduction of five-dimensional i?^-Lagrangians to four 



dimensions [3|-t5|. Similar non- minimal coupling terms were obtained by calculating QED 
one-loop polarization on a curved background (6|. Recently, the behavior of the rotational 
velocities of test particles gravitating around galaxies was investigated in modified gravity 
for non-minimal matter couplings jr-Q]. Furthermore, the non-minimal couplings of gravity 
and electromagnetism were extended to I{R)F^ form to explain late-time cosmic acceleration 
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llL primordial magnetic field during the reheating epoch 12| and rotation curves of 



galaxies [13 1 . 



Nevertheless, the non-minimal couplings have not been investigated in sufficient detail 
as regards to applications such as the contributions to the rotational curve of galaxies and 
Pioneer anomaly. To gain more insights on such configurations in non-minimal models, we 
proceed to investigate exact solutions of them for spherically symmetric systems with charge 
and mass. In order to find more general solutions which are compatible with observations 
from the solar system to cosmological scales, we propose i?''F^-type couplings between 
gravity and electromagnetism through arbitrary real power of the curvature scalar. We 
determine the field equations by a first order variational principle using the method of 
Lagrange multipliers. Then we present a wide range of exact static, spherically symmetric 
solutions that depend on various values of (3. It is interesting to note that the solutions are 
asymptotically fiat in general for values of /3 7^ within the range — | < /3 < 1. 

The non-minimal couplings in general give contributions both to the Maxwell and the 
Einstein field equations. Contributions to the Maxwell equations can be associated with 
the magnetization and the polarization of a specific medium while contributions to the 
Einstein equations may give important modifications to the space-time metric. Such effects. 



if there are any, may be used to explain some unexpected observations of gravity such as 
dark matter, dark energy and Pioneer anomaly considering only the electromagnetic and 
gravitational fields. 



II. FIELD EQUATIONS OF THE EXTENDED EINSTEIN-MAXWELL THEORY 



We will obtain our field equations by a variational principle from an action 

J[e^a;^,F]= / L = [ C*l (1) 
Jm Jm 

where {e"} and {uj""^} are the fundamental gravitational field variables and F is the electro- 
magnetic field 2-form. The space-time metric g = r/a^e" ® with signature ( — h ++) and 
we fix the orientation by setting *1 = e° A A A e^. Torsion 2-forms and curvature 
2-forms i?"^ of space-time are found from the Cartan-Maurer structure equations 

de"" + uj\ A e'' = T% (2) 

duj\ + u\ A u\ = R\. (3) 

We consider the following Lagrangian density 4-form which involves R^F'^ terms to all orders 
in /3: 

1 1 °° 

L = —R*l--Y.{{a,RfTFA*F (4) 

where is Newton's universal gravitational coupling constant, ao is a coupling constant with 

dimension [L]^ and /3 is a real number. When we take ao = 0, we get back to the minimal 

Einstein-Maxwell theory. If we take ao 7^ and assume that the condition \{aQRY\ < 1 is 

00 

satisfied, we can write {{o-oR) = i-(aoR)f^ • '^^^ right hand side diverges as /3 — t- or 



n=0 



— )■ ^. Therefore, the values (3 = and ao-R = 1 should be avoided. Thus we replace the 



Lagrangian (jl]) with 



L = ^R*l - — 1 -^FA*F. (5) 

2«;2 2(1 - {aoRY) ^ ' 



which is meaningful even when the condition |(ao-R)^| < 1 is not satisfied. 

The field equations are obtained by considering the independent variations of the action 
with respect to {e"}, {<^'*^} and {F}. The electromagnetic field components are read from the 



expansion F — \Fabe°' A e^. We will choose the unique metric-compatible, torsion-free Levi- 
Civita connection. We impose this choice of the connection through constrained variations 
by the method of Lagrange multipliers. We constrain the electromagnetic field 2-form F to 

be closed, that is; dF = 0, and this is imposed by the variation of a Lagrange multiplier 
2-form /i. Therefore, we add to the above Lagrangian density the following constraint terms: 

Lc = (de"" + A e^) A A„ + A (6) 

where Aa's are Lagrange multiplier 2- forms whose variation imposes the zero-torsion con- 
straint T" = 0. 

The infinitesimal variations of the total Lagrangian density L -(- Lc (modulo a closed form) 
is found to be 

L + Lc^ Are" A R^^ A *e„bc + e" A — 1 r^i^aF A *F - F A t„ * F) + e" A DX^ 

Ba^B^~^ 1 
+e" A ^i_l^^jiyy ('aR')(ibF A*F + FAl,*F) + -Cj^ A (e^ A A" - e" A A^) 

+cO„6AE«^-FA^— ^^-^*F + AaAT"-FAd^ (7) 

where the angular momentum tensor E"'': 



1 fia^Rf^- 
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^ — F)f-— ^ -^{F"^ * F + F^ At"" * F - F'' Ai^ * F - F Ai""^ * F)]. (8) 

2 ^(1 - {aoRyy^ ^ ' 

We will use the abbreviations A A • • • = e"**'", UaF — Fa, ihaF — Fab, i^aR^'b — 
Rb, LbaR""^ — R- Lagrange multiplier 2-forms A^ are solved from the connection variation 
equations 

Ca A Ab - 66 A A„ = 2Y.ab (9) 

by applying the interior product operators twice as 

A" = 2^bE^'^ -I- \ibc^'^ A e\ (10) 

When we substitute the Aa's above into the co-frame equation, we find the Einstein field 
equations for the non-minimal theory: 

^^i?"'^ A *e„6e + ^^^-J— ^(,,F A *F - F A i„ * F) + D{2ibY.'>- + ^e" A ^^.eE'^^) 



^a^R^-' 

{l-{aoRy 



^{iaR'){ibFA*F + FAib*F)^0. (11) 



The Maxwell equations read 

We may represent the effects of non-minimal couplings of the electromagnetic fields to gravity 
through the definition of a constitutive tensor. Thus, Maxwell's equations for an electro- 
magnetic field F in an arbitrary medium can be written as 

dF = {) , *d*G=J (13) 

where G is called the excitation 2-form and J is the source electric current density l-form. 
The effects of gravitation and electromagnetism on matter are described by G and J. We 
can complete this system using electromagnetic constitutive relations relating G and J to 
F. Here we consider only the source-free interactions, that is J = 0. Then, we can write a 
simple linear constitutive relation 

G = Z{F) (14) 
where ^ is a type-(2,2)-constitutive tensor. For the above theory, we have 

We can identify the polarization l-form p — i^(aQR)i^ ^u-^ and the magnetization l-form 
m = — xl"(a^fl)/3 ''U * F related with the above theory where U is a unit, time-like velocity 
vector field associated with an inertial observer. One can find more information about these 



concepts in 14\- 



III. STATIC, SPHERICALLY SYMMETRIC SOLUTIONS 

We look for static, spherically symmetric solutions to the field equations which are given by 
the metric 

g = -f{rfde + f{r)-^dr^ + r^dO^ + sin2(6')#2. (i6) 
We consider a static electric potential l-form A = h{r)dt. Then, electromagnetic field 2-form 

F = dA = h'dr A dt = Hdr A dt (17) 
has a spherically symmetric electric field component only. 



We note that, for the minimal case ao = 0, the Reissner-Nordstrom metric 

^ = _(1 _ ^ + 'f^)df + (1 - ^ + ^)-Mr2 + r'de' + sm'{e)dcP^ (18) 
and the Coulomb electric potential 

A = -'^dt (19) 
r 

is a solution of f lTTj) and f|T2|) . 

After a lengthy calculation, the non-minimally coupled (ao 7^ 0) Einstein-Maxwell equations 
ffTTj) and f|T2|) are reduced for the metric f|T6l) and the electromagnetic 2-form f|T7|) . We do 
not find it useful to write down the reduced field equations in full detail, because they are 
very long and complicated. In order to avoid the problems of having higher order derivatives 
in our field equations, we restrict attention to those cases for which the Lagrange multipliers 
Aq = 0. We calculate 

3. y^,«=_M5!__), (20) 



where the curvature scalar is given by 



-i? = f" + ^ + |(f -1). (21) 



Thus, we set 

The Maxwell equations (fT2|) can be immediately integrated to give 

1 



= constant. (22) 



1 - {aoRf 

where q is the electric charge determined by the Gauss integral 



Hr' = q (23) 



*G = q. 



The expressions fl22l) and fl23l) imply the following relation for consistency: 

R = a^r A (24) 

where ai is an integration constant to be fixed. Then, it can be shown easily that the 
remaining Einstein field equations become 



R,^ K^(3a^oR^'^ . 



The coupled equations (12^ -( 126|) are solved by 

rlr . C ^2^2 a (/3 _ 1)2 2,3+2 ^ , n , 1 

, / N 9 — 1) 3/3 + 1 

with tti = («;2g2/3ao)~. go back to non-minimally coupled Einstein-Maxwell field 

equations and check the solutions, we fix the integration constant in the fl22|) to be We 
further note that the values of /3 7^ in the range — | < /3 < 1 give asymptotically flat 
solutions. For /? in this interval, we identify the mass M = Moreover, r = is an 
essential singularity in general where the quadratic curvature invariant *{Rab A *R"-^) — )■ 00 
as r — )■ 0. It is remarkable to observe that the solution (1271) goes to the Schwarzschild 
solution for q = 0, while for q but uq = the metric in fl27j) does not coincide with the 
Reissner- Nordstrom metric (fT8|) . The explanation is as follows: For uq = 0, we have R = 0, 
so that ( l26i) is identically satisfied. In this case, (1251) is solved by the Reissner- Nordstrom 
metric. On the other hand, for ao 7^ and ai 7^ we solve (l26l) with the choice 



(l-(aoW 

for which the limit ao — > does not exist. We would like to point out that with the choice 



H^ = 



[I - (aoR) 



/3^2 



for an arbitrary ag 7^ 0, we may set ai = in and show that our non-minimally coupled 
model admits the Reissner-Nordstrom black hole (fT8|) and (fT9l) as a solution. 

A metric function similar to ours above with /3 = | is obtained and analyzed in a recent 
paper 15|. They consider a different theory of non-minimally coupled electromagnetic fields 
to gravity. Furthermore, their electromagnetic field is that of a Dirac magnetic monopole 
Qm- Contrary to our case, this solution joins smoothly to the Reissner-Nordstrom solution 
in the absence of non-minimal couplings. We also wish to point out that the metric fl7r|) with 
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3 



in a dilaton-gravity 



(3 = —3 involves a Rindler acceleration term. This is noted by 
model and may be related with the source of various anomalies such as the rotation curves 
of spiral galaxies and the Pioneer anomaly. 



IV. CONCLUSION 



We have formulated non-minimal i?^F^-type coupled Einstein-Maxwell theory. The field 
equations are obtained by a first order variational principle using the method of Lagrange 
multipliers in the language of exterior differential forms. We found static, spherically sym- 
metric solutions for arbitrary values of the (3 parameter. Such solutions may be related 
with the source of various anomalies in gravity for certain values of the parameters Uq 
and /3. Moreover, they may contribute to resolution of many important challenges such as 
dark matter, dark energy and Pioneer anomaly without any other fields. In addition, the 
non-minimal terms modify electromagnetic potentials at cosmological scales. Then the con- 
ventional electromagnetic energy density of the universe gets modified. In other words; if the 
effects of dark matter are not due to some exotic matter fields, the non-minimal couplings 
can lead to such effects[18:]. Even if the electric charge q is not large, one may choose the 
parameters in such a way that these solutions can explain the rotation curves of galaxies 
and Pioneer anomaly for some parameter values. 
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